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Abstract 

Angular correlations between unidentified charged trigger (t) and associated (a) particles 
are measured by the ALICE experiment in Pb-Pb collisions at -^/^nn = 2.76 TeV for trans- 
verse momenta 0.25 < pj-" < 15 GeV/c, where pj > pf. The shapes of the pair correlation 
distributions are studied in a variety of collision centrality classes between and 50% of 
the total hadronic cross section for particles in the pseudorapidity interval | T] | < 1.0. Dis- 
tributions in relative azimuth A0 = 0' — (p" are analyzed for |At]| = Itj' — T]"| > 0.8, and 
are referred to as "long-range correlations". Fourier components V„a = (cos(?iA0)) are ex- 
tracted from the long-range azimuthal correlation functions.If particle pairs are correlated to 
one another through their individual correlation to a common symmetry plane, then the pair 
anisotropy VnA{pj,pf) is fully described in terms of single-particle anisotropics v„{pj) as 
yn\{PT^PT) ~ Vn{PT)vn{PT)- This expectation is tested for 1 < « < 5 by applying a global 
fit of all V„A{Px,pf) to obtain the best values v„{GF}(/7r). It is found that for 2 < « < 5, 
the fit agrees well with data up to pj ~ 3-4 GeV/c, with a trend of increasing deviation as 
p'j. and pf are increased or as collisions become more peripheral. This suggests that no 
pair correlation harmonic can be described over the full 0.25 < pr < 15 GeV/c range using 
a single Vnipr) curve; such a description is however approximately possible for 2 < « < 5 
when pj <4 GeV/c.For the n = l harmonic, however, a single vi (pr) curve is not obtained 
even within the reduced range pf <4 GeV/c. 



*See Appendix A for the list of collaboration members 
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1 Introduction 

Ultra-relativistic collisions of large nuclei at the Large Hadron Collider (LHC) and at the Rel- 
ativistic Heavy Ion Collider (RHIC) enable the study of strongly-interacting nuclear matter at 
extreme temperatures and energy densities. One key piece of evidence for the formation of 
dense partonic matter in these collisions is the observation of particle momentum anisotropy in 
directions transverse to the beam [[T]-[6l. One powerful technique to characterize the properties 
of the medium is with two-particle correlations [|7]-fT8l. which measure the distributions of an- 
gles and/or Arj between particle pairs consisting of a "trigger" at transverse momentum p'j. 
and an "associated" partner at pf. 

In proton-proton collisions, the full (A^, At]) correlation structure at (A(^, At]) ^ (O5O) is 
dominated by the "near-side" jet peak, where trigger and associated particles originate from 
a fragmenting parton, and at A^ ^ :;r by the recoil or "away-side" jet. The away-side peak is 
broader in At] , due to the longitudinal momentum distribution of partons in the colliding nu- 
clei. In central nucleus-nucleus collisions at RHIC, an additional "ridge" feature is observed 
at A^ ^ fl 13111411 . which has generated considerable theoretical interest [fT9l - |29ll since its ini- 
tial observation. With increasing px, the contribution from the near-side jet peak increases, 
while the ridge correlation maintains approximately the same amplitude. The recoil jet correla- 
tion is significantly weaker than that of the near side, because of kinematic considerations [|30ll 
and also because of partonic energy loss. When both particles are at high transverse momenta 
(Pt ~ 6 GeV/c), the peak shapes appear similar to the proton-proton case, albeit with a more 
suppressed away side. This away-side correlation structure becomes broader and flatter than in 
proton-proton collisions as the particle pr is decreased. In fact, in very central events (~ 0- 
2%), the away side exhibits a concave, doubly-peaked feature at |A^ — :;r| ~ 60° [|3T| , which 
also extends over a large range in |At] | [fTTl [TSl . The latter feature has been observed previ- 
ously at RHIC [fT2l - [l4l . but only after subtraction of a correlated component whose shape was 
exclusively attributed to elliptic flow. 

However, recent studies suggest that fluctuations in the initial state geometry can generate 
higher-order flow components |32]-|40ll. The azimuthal momentum distribution of the emitted 
particles is commonly expressed as 

dN 

— ocl + £ 2Vn{pT) COS (1) 
' n=l 

where v„ is the magnitude of the order harmonic term relative to the angle of the initial- state 
spatial plane of symmetry First measurements, in particular of V3 and V5 have been reported 
recently El [SB El. 

These higher-order harmonics contribute to the previously-described structures observed in 
trigger-associated particle correlations via the expression 

o^l + Y, 2v„(/r)v„br) COS (nA0) . (2) 

Similarly, the measured anisotropy from two-particle correlations at harmonic order n is given 
by V„a: 

HNPairs 

oc 1 + £ 2V„A(/r,Pr) COS {nA^) . (3) 
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In this article, we present a measurement of the V„a coefficients from triggered, pseudorapidity- 
separated (|A77| > 0.8) pair azimuthal correlations in Pb-Pb collisions in different centrality 
classes and in several transverse momentum intervals. Details of the experimental setup and 
analysis are described in sections |2] and [3| respectively. The goal of the analysis is to quanti- 
tatively study the connection between the measured two-particle anisotropy VnA of Eq. |3] and 
the inclusive-particle harmonics of Eq. [2j Specifically, we check whether a set of single- valued 
Vnipr) points can be identified that describe the measured long-range anisotropy via the relation 
Vn{pj)vn{pj) = V„/^{pj,pj). If so, VnA is Said to factorize into single-particle Fourier coeffi- 
cients within the relevant Pj^Px region. This relationship is tested for different harmonics n 
and in different centrality classes by performing a global fit (GF) over all pj" bins (see sec- 
tion |4]). The global fit procedure results in the coefficients v„{GF}(pr) that best describe the 
anisotropy given by the VnAip'j-^Pj) harmonics as Vn{GF}{pj) x v„{GF}{pj). The resulting 
Vn{GF} values for 1 < n < 5 are presented in section [5] A summary is given in section |6] 

2 Experimental setup and data analysis 

The data used in this analysis were collected with the ALICE detector in the first Pb-Pb run 
at the LHC (November 2010). Charged particles are tracked using the Time Projection Cham- 
ber (TPC), whose acceptance enables particle reconstruction within — 1.0 < 7] < 1.0. Primary 
vertex information is provided by both the TPC and the silicon pixel detector (SPD), which con- 
sists of two cylindrical layers of hybrid silicon pixel assemblies covering | rj | < 2.0 and 1 7] | < 1 .4 
for the inner and outer layers, respectively. Two VZERO counters, each containing two arrays of 
32 scintillator tiles and covering 2.8 < 7] < 5.1 (VZERO-A) and -3.7 < 7] < -1.7 (VZERO- 
C), provide amplitude and time information for triggering and centrality determination. The 
trigger was configured for high efficiency to accept inelastic hadronic collisions. The trigger is 
defined by a coincidence of the following three conditions: i) two pixel hits in the outer layer of 
the SPD, ii) a hit in VZERO-A, and iii) a hit in VZERO-C. 

Electromagnetically induced interactions are rejected by requiring an energy deposition above 
500 GeV in each of the Zero Degree Calorimeters (ZDCs) positioned at ± 114 m from the 
interaction point. Beam background events are removed using the VZERO and ZDC timing 
information. The combined trigger and selection efficiency is estimated from a variety of Monte 
Carlo (MC) studies. This efficiency ranges from 97% to 99% and has a purity of 100% in the 
0-90% centrality range. The dataset for this analysis includes approximately 13 million events. 
Centrality was determined by the procedure described in Ref . [42J . The centrality resolution, 
obtained by correlating the centrality estimates of the VZERO, SPD and TPC detectors, is 
found to be about 0.5% RMS for the 0-10% most central collisions, allowing centrality binning 
in widths of 1 or 2 percentiles in this range. 

This analysis uses charged particle tracks from the ALICE TPC having transverse momenta 
from 0.25 to 15 GeV/c. The momentum resolution o{pt) / Pt rises with px and ranges from 1- 
2% below 2 GeV/c up to 10-15% near 15 GeV/c, with a negligible dependence on occupancy. 
Collision vertices are determined using both the TPC and SPD. Collisions at a longitudinal po- 
sition greater than 10 cm from the nominal interaction point are rejected. The closest-approach 
distance between each track and the primary vertex is required to be within 3.2 (2.4) cm in the 
longitudinal (radial) direction. At least 70 TPC pad rows must be traversed by each track, out 
of which 50 TPC clusters must be assigned. In addition, a track fit is applied requiring per 
TPC cluster < 4 (with 2 degrees of freedom per cluster). 
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Fig. 1: Examples of two-particle correlation functions C(A0, Atj) for central Pb-Pb collisions at low to 
intermediate transverse momentum (left) and at higher pj (right). Note the large difference in vertical 
scale between panels. 

3 Two-particle correlation function and Fourier analysis 

The two-particle correlation observable measured here is the correlation function C(A^, At]), 
where the pair angles A^ and At] are measured with respect to the trigger particle. The correla- 
tions induced by imperfections in detector acceptance and efficiency are removed via division 
by a mixed-event pair distribution Admixed (^^ 5 in which a trigger particle from a particular 
event is paired with associated particles from separate events. This acceptance correction pro- 
cedure removes structure in the angular distribution that arises from non-uniform acceptance 
and efficiency, so that only physical correlations remain. Within a given pj, pj, and centrality 
interval, the correlation function is defined as 

C(A<^,AT])^^x ^— (4) 

^ " A^same Admixed ( A(^), AT]) 

The ratio of mixed-event to same-event pair counts is included as a normalization prefactor 
such that a completely uncorrected pair sample lies at unity for all angles . For Admixed i^^i^l)^ 
events are combined within similar categories of collision vertex position so that the acceptance 
shape is closely reproduced, and within similar centrality classes to minimize effects of residual 
multiplicity correlations. To optimize mixing accuracy on the one hand and statistical precision 
on the other, the event mixing bins vary in width from 1 to 10% in centrality and 2 to 4 cm in 
longitudinal vertex position. 

It is instructive to consider the two examples of C(A^, At]) from Fig. [T] to be representative of 
distinct kinematic categories. The first is the "bulk-dominated" regime, where hydrodynamic 
modeling has been demonstrated to give a good description of the data from heavy-ion colli- 
sions [IIH3. We designate particles with pj (thus also pj) below 3-4 GeV/c as belonging to 
this region for clarity of discussion (see Fig.[T] left). A second category is the "jet-dominated" 
regime, where both particles are at high momenta (p^ > 6 GeV/c), and pairs from the same 
di-jet dominate the correlation structures (see Fig.[T| right). 

A major goal of this analysis is to quantitatively study the evolution of the correlation shapes 
between these two regimes as a function of centrality and transverse momentum. In order to 
reduce contributions from the near-side peak, we focus on the correlation features at long range 
in relative pseudorapidity by requiring |At] | > 0.8. This gap is selected to be as large as possible 
while still allowing good statistical precision within the TPC acceptance. The projection of 
C(A<^, |At] I > 0.8) into A(j) is denoted as C(A(^). 
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Pb-Pb 2.76 TeV, 0-2% central 
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Fig. 2: (Color online) Left: C (A0) for particle pairs at | At] | > 0.8. The Fourier harmonics for Via to Vsa 
are superimposed in color. Their sum is shown as the dashed curve. The ratio of data to the n <5 sum is 
shown in the lower panel. Center: Amplitude of VnA harmonics vs. n for the same pj, pj, and centrality 
class. Right: V„a spectra for a variety of centrality classes. Systematic uncertainties are represented with 
boxes (see section]?]), and statistical uncertainties are shown as error bars. 
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Fig. 3: (Color online) Left: C(A(^) at |At]| > 0.8 for higher-pr particles than in Fig. ]2] The Fourier 
harmonics VnA for « < 5 are superimposed in color. Their sum is shown as the dashed curve. The ratio of 
data to the n<5 sum is shown in the lower panel. Right: Amplitude of V„a harmonics vs. n at the same 
Pj, pj for two centrality bins. Systematic uncertainties are represented with boxes (see section ]4]l, and 
statistical uncertainties are shown as error bars. 



An example of C(A(^) from central Pb-Pb collisions in the bulk-dominated regime is shown 
in Fig. |2] (left). The prominent near-side peak is an azimuthal projection of the ridge seen in 
Fig. JT] In this very central collision class (0-2%), a distinct doubly-peaked structure is visible 
on the away side, which becomes a progressively narrower single peak in less central colli- 
sions. We emphasize that no subtraction was performed on C(A^), unlike other jet correlation 
analyses HTHT?!. 

A comparison between the left panels of Fig. |2]and Fig. ]3] demonstrates the change in shape 
as the transverse momentum is increased. A single recoil jet peak at ~ ;r appears whose 
amplitude is no longer a few percent, but now a factor of 2 above unity. No significant near-side 
ridge is distinguishable at this scale. The recoil jet peak persists even with the introduction of a 
gap in I At] I due to the distribution of longitudinal parton momenta in the colliding nuclei. 

The features of these correlations can be parametrized at various momenta and centralities by 
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decomposition into discrete Fourier harmonics, as done (for example) in BSl I40B . Following 
the convention of those references, we denote the two-particle Fourier coefficients as (see 
Eq. [3]), which we calculate directly from C (A<^) as 



V„A = (cos (nA0)) = ^C/COs(nA0/) / . 



(5) 



Here, Ci indicates that the C (A0) is evaluated at A^/. Thus y„A is independent of the normaliza- 
tion of C (A0). The VnA harmonics are superimposed on the left panels of Fig. [2] and Fig. |3j In 
the right panels, the V^a spectrum is shown for the same centrality and momenta, with additional 
centrality classes included to illustrate the centrality dependence. The systematic uncertainties 
in these figures are explained in section |4j 
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Fig. 4: V„A coefficients as a function of /j'^ for the 0-2%, 0-10%, and 40-50% most central Pb-Pb 
collisions (top to bottom). 

In the bulk-dominated momentum regime and for central collisions (Fig. [2]), the first few Fourier 
harmonics are comparable in amplitude, with the notable exception of Via- The first 5 combined 
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harmonics reproduce C(A0) with high accuracy, as shown in the ratio between the points and 
the component sum. For less central collisions, V2A increasingly dominates. In the high-p^ 
regime (Fig. [3]), the jet peak at = ;r is the only prominent feature of the correlation function. 
The even (odd) harmonics take positive (negative) values which diminish in magnitude with 
increasing n, forming a pattern distinct from the low-pr case. The dependence of the values 
on n in the left panel of Fig. |3]is approximately consistent with a Gaussian function centered 
at n = 0, as expected for the Fourier transform of a Gaussian distribution of width C7„ = 1 / Oa(^ 
centered at A^ = n. In this case, the sum of the first 5 harmonics does not reproduce C(A^) 
with the accuracy of the low-pj- case, as suggested by the larger value (61.5/35 compared to 
33.3/35). Although not shown, it was found that including additional harmonics significantly 
improves the X^/NDF measure for high-p^ correlations, but adding higher orders in the bulk- 
dominated case has only a modest effect. For example, if curves composed of the lowest ten 
harmonics are used, the value drops by only about 10% to 30.0/35 in Fig. [2} but by over 
25% to 45.7/35 in Fig. |3j We note that V2 is not the dominant coefficient in Fig. [3| instead, its 
magnitude fits into a pattern without significant dependence on collision geometry, as suggested 
by the continuous decrease with increasing n for both 0-20% and 40-50% central events. This 
suggests that the n spectrum is driven predominantly by intra-jet correlations on the recoil side, 
as expected from proton-proton correlations at similar particle momenta. 

Figure |4] shows the V„a coefficients as a function of trigger pr for a selection of associated pr 
values. For n >2, VnA reaches a maximum value at pj ~ 3-4 GeV/c, decreasing toward zero 
(or even below zero for odd n) as p'j. increases. This rapid drop of the odd coefficients at high 
provides a complementary picture to the n dependence of V^a shown in Fig. [3] 

4 Factorization and the global fit 

The trends in pj and centrality in Fig. |4] are reminiscent of previous measurements of v„ from 
anisotropic flow analyses [[171 [3T1 [411. This is expected if the azimuthal anisotropy of final 
state particles at large jAr/ 1 is induced by a collective response to initial-state coordinate-space 
anisotropy from collision geometry and fluctuations [[37l . In such a case, C(A(^) reflects a 
mechanism that affects all particles in the event, and VnA depends only on the single-particle 
azimuthal distribution with respect to the n-th order symmetry plane Under these circum- 
stances VnA factorizes as 

= (v„{2}(pV)v„{2}(p^)). (6) 

Here, () indicates an averaging over events, (()) denotes averaging over both particles and 
events, and v„{2} specifies the use of a two-particle measurement to obtain v„. 

Equation[6]represents the factorization ofV„A{p'f,Pj) into the event-averaged product {v„{2}{pj) 
Vn{2}{pj)), which includes event-by-event fluctuations. Consistency with Eq. |6] suggests that 
a large fraction of the particle pairs are correlated through their individual correlation with a 
common plane of symmetry. For example, symmetry planes for particle pairs at pj and pj- may 
develop at a harmonic order n from collision geometry, initial state density fluctuations, or from 
an axis formed by (di-)jet fragmentation. If a single- valued Vnipr) curve on an interval con- 
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taining p^j and pj can reproduce the magnitude of any Vnh{pj,Pj), then y„A factorizes within 
the (pj.Pj) region. 

Equation [6] is tested by applying a global fit to the V„a data points over all pj and pj bins 
simultaneously. This is done separately at each order in n and for each centrality class. An 
example from the 0-10% most central event class is shown in Fig. [5| where the VnA points for 
n = 2 to 5 are plotted (in separate panels) on a single pj, pj axis as indicated. The global 
fit function depends on a set of unconstrained and independent parameters, where N is the 
number of pj (or pj) bins. The parameters are v„{GF}(pr), with the fit generating the product 
Vn{GF}{pj) X Vn{GF}{pj) that minimizes the total for Ka points. 

The sources of systematic uncertainty of Vn^ are those that cause A0-dependent variation on 
C(A^). Factors affecting overall yields such as single-particle inefficiency cancel in the ratio 
of Eq. |4| and do not generate uncertainty in C(A^). Table [T] shows the different contributions 
to the systematic uncertainty of Vn^, and Table |2] lists typical magnitudes of these uncertainties 
for a few representative centrality classes. 

The event mixing uncertainty (denoted as "a" in Table [T]) accounts for biases due to imperfect 
matching of event multiplicity and collision vertex position, as well as for finite mixed-event 
statistics. This uncertainty changes with pj, pj, and centrality. It is evaluated by comparing the 
n<5 Fourier sum from C(A0) with that from A^same(A^). The uncertainty from (a) is depicted 
by grey bars on the points in C (A^) in Fig. [2] and Fig.|3j Due to fluctuations in the mixed-event 
distribution, uncertainty (a) tends to scale with the y„A statistical error, as shown in the table. 

The remainder of the systematic uncertainties are not assigned to C (A0), but rather to each y„A 
directly, where their influence on y„A is more clearly defined. The uncertainty from centrality 
determination (b) accounts for the resolution and efficiency of the detector used for multiplicity 
measurements, as well as any biases related to its r\ acceptance. This uncertainty is globally 
correlated in centrality. It was studied by conducting the full analysis with the SPD as an alter- 
native centrality estimator, since it has different systematic uncertainties and covers a different 
pseudorapidity range than the VZERO detectors. The results were found to agree within 1%. 

Uncertainty from tracking and momentum resolution (c) was evaluated on C (A0) using differ- 
ent track selection criteria. Slightly larger correlation strength is obtained for more restrictive 
track selection (at the expense of statistical loss), and the difference was found to grow with pj 
by roughly 1% per GeV/c. 

Additional uncertainty is introduced by (d) the finite A0 bin width, which was estimated by 
comparing the RMS bin width to the harmonic scale, and (e) the precision of the extraction 
(Eq. [5]). The latter was estimated by calculating (sin(nA(^)), which is independent of n, and 
should vanish by symmetry. The residual finite values are used to gauge the corresponding 
y„A uncertainty. Because the amplitude of the y„A harmonics tends to diminish with increasing 
n, both uncertainties (d) and (e) are small for n < 6 but become comparable to y„A at higher n. 
Effects (a)-(e) are all combined in quadrature to produce the y„A systematic uncertainties, which 
are depicted as the solid colored bars on the points in Fig.|4]and Fig.|5j Finally, the uncertainty 
(f) is included in the quadrature sum with (a)-(e) for yiA only. 
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Fig. 5: (Color online) Global fit examples in 0-10% central events for n = 2,3,4 and 5. The measured 
VnA coefficients are plotted on an interleaved pj, pj axis in the upper panels, and the global fit function 
(Eq. |6]l is shown as the red curves. The global fit systematic uncertainty is represented by dashed lines. 
The lower section of each panel shows the ratio of the data to the fit, and the shaded bands represent the 
systematic uncertainty propagated to the ratio. In all cases, off-scale points are indicated with arrows. 
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Contribution Magnitude 

(a) Event mixing 20-30% astat 

(b) Centrality determination 1% 

(c) Track selection, pj resolution l%y„A x {Pj) 

(d) bin width (0.8n)%y„A 

(e) y„A extraction < 10% Vn^ (n<6); 10-30% V'„a(«>6) 

(f) Vn^ {n=\ only) 10%y„A x 

Table 1: Systematic uncertainties on Vn\- 
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Centrality 


(Ka) 


{^sys{tot)l 




Ob 






Oe 








(xlO-3) 


(Xl0~3) 






(xlO-3) 








1 


0-10% 


4.2 


13 


0.57 


0.042 


0.042 


0.033 


13 


0.42 


1 


20-30% 


11 


3.7 


0.6 


0.11 


0.11 


0.089 


3.5 


1.1 
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40-50% 
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0.38 
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0-10% 


12 
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0.19 
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31 
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0.46 


0.31 


0.31 


0.5 


0.6 





2 


40-50% 


43 


6.9 


0.4 


0.43 


0.43 


0.69 


6.8 





3 


0-10% 


0.82 


18 


0.29 


0.0082 


0.0082 


0.02 


18 





3 


20-30% 


2.1 


1.4 


0.42 


0.021 


0.021 


0.05 


1.4 





3 


40-50% 


10 


4.2 


0.42 


0.1 


0.1 


0.25 


4.2 





4 


0-10% 


2.6 


5.4 


0.44 


0.026 


0.026 


0.083 


5.4 





4 


20-30% 


7.2 


0.45 


0.27 


0.072 


0.072 


0.23 


0.26 





4 


40-50% 


11 


2.6 


0.47 


0.11 


0.11 


0.35 


2.6 





5 


0-10% 


2.6 


6.2 


0.32 


0.026 


0.026 


0.1 


6.1 





5 


20-30% 


2.8 


3.2 


0.17 


0.028 


0.028 


0.11 


3.2 





5 


40-50% 


6.8 


0.6 


0.28 


0.068 


0.068 


0.27 


0.45 






Table 2: Typical values of Vn^ systematic uncertainties. 



To evaluate the systematic uncertainty, the global fit procedure is performed three times for 
each n and centrality bin: once on the measured y„A points (leading to the red curves in Fig.|5]), 
and once on the upper and lower bounds of the systematic error bars (resulting in black dashed 
curves). The v„{GF} systematic error is then assigned as half the difference. The resulting 
uncertainties are shown as open boxes in Fig.|6]and Fig. 1 1 , which are discussed in the following 
sections. 



5 Global fit results 

In the n = 2 case (Fig. [sj top), the fit agrees well with the data points at low and pj, 
but diverges with increasing pj for each pj interval. Where disagreement occurs, the fit is 
systematically lower than the points. In contrast, for n = 3, the fit does not follow the points 
that drop sharply to negative values at the highest momenta. This is also observed for n = 5, 
though with poorer statistical precision. 

The global fit is driven primarily by lower particle pt, where the smaller statistical uncertainties 
provide a stronger constraint for minimization. The disagreement between data and the fit, 
where pj and pj are both large, points to the breakdown of the factorization hypothesis; see 
also Fig. [3] and the accompanying discussion. 
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The factorization hypothesis appears to hold for n > 2 at low pj (< 2 GeV/c) even for the 
highest p'j. bins. The y„A values for these cases are small relative to those measured at higher 
Pj, and remain constant or even decrease in magnitude as p'j- is increased above 3-4 GeV/c. 
V2A dominates over the other coefficients, and the n> 3 terms are not significantly greater than 
zero. This stands in contrast to the high-p^, high-py case, where it was demonstrated in Fig.|3] 
that dijet correlations require significant high-order Fourier harmonics to describe the narrow 
recoil jet peak. 

0.25 
0.2 
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O 0.1 
"Z=0.05 
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« • o o • 

8 o o o 
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0.25<p^<15GeV 



.1'" ' f • 
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|S i 



f 



HI 

1- 



[GeV/c] 

Fig. 6: (Color online) The global-fit parameters, v„{GF}, for 2 < « < 5. Statistical uncertainties are 
represented by error bars on the points, while systematic uncertainty is depicted by open rectangles. 

The parameters of the global fit are the best- fit Vn{GF} values as a function of pj, which can 
be interpreted as the coefficients of Eq. [T] The results of the global fit for 2 < n <5, denoted 
Vn{GF}, are shown in Fig. [6] for several centrality selections. We note that the global fit con- 
verges to either positive or negative v„{GF} parameters, depending on the starting point of the 
fitting routine. The two solutions are equal in magnitude and goodness-of-fit. The positive 
curves are chosen by convention as shown in Fig. [6} In the 0-2% most central data, vsjGF} 
(v4{GF}) rises with pj relative to V2{GF} and in fact becomes larger than V2{GF} at ap- 
proximately 1.5 (2.5) GeV/c. V2{GF} reaches a maximum value near 2.5 GeV/c, whereas the 
higher harmonics peak at higher pr- These data are in good agreement with recent two-particle 
anisotropic flow measurements \3U at the same collision energy, which included a pseudora- 
pidity gap of |At] | > 1.0. 

For 2 < n < 5, the results are not strongly sensitive to the upper p^ limit included in the global 
fit. The global fit was performed not only over the full momentum range (as shown in Fig. |6]), 
but also with the restriction to Vha points with pf < 2.5 GeV/c. The outcome was found to 
be identical to the full fit within one standard deviation. This again reflects the weighting by 
the steeply-falling particle momentum distribution, indicating that a relatively small number of 
energetic particles does not strongly bias the event anisotropy, as calculated by the global fit. 

If the global fit is applied to y„A points exclusively at large particle momenta, factorization 
behavior can be tested for correlations that are predominantly jet-induced. An example is shown 
in figure Fig.|7| where the global fit has been applied to VnA points within 5 < pf < 15 GeV/c. In 
this case, there are six VnA datapoints fitted, and tliree fit parameters, which are v„{GF} at 5-6, 
6-8, and 8-15 GeV/c. An approximate factorization is observed over this range. The agreement 
between fit and data for the lowest fitted datapoint (at 5-6 GeV/c) is rather poor, indicating that 
the correlations there are in a transitional region that is less jet-dominated than at higher px- 

The parameters from the high-pj- global fit can be plotted, just as was done in Fig.[6[ to demon- 
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Fig. 7: (Color online) High-^^r fit examples in 0-20% central events for « = 1 to 4. Although all 
datapoints are shown for p'j, > 5 GeV/c, the fit range includes only the six points with Pj>5 GeV/c. 

strate their pj and centrality dependence. However, the sign definition of v„ becomes problem- 
atic in the case where v„{pj) and Vn{pj) have the same sign, but y„A < 0- In this case, the v„ 
coefficients are represented to be positive as a matter of convention. 

The fit results from these high-p^, high-pf- long-range correlations are shown as open points in 
Fig. [8] The clear deviation between the two different sets of points demonstrates that it is not 
possible for a single-valued set of v„(pr) points to simultaneously describe both low-p^ and 
high-p^ pair anisotropy. 

1 

0.8 

_ 0.6 

i 0.4 
> 

0.2 


-0.2 

U O lU U O lU u 

[GeV/c] 

Fig. 8: (Color online) The global-fit parameters, v„{GF}, for 1 < « < 5 as obtained using restricted p^ 
fit ranges at two different centralities. The solid (open) points represent fits using only 0.25 < pj < 5 
(5 < pj < 15) GeV/c. The open points represent the magnitude of v„{GF} from high-/7y, high-/7y long- 
range correlations. Statistical uncertainties are represented by error bars on the points, while systematic 
uncertainty is depicted by open rectangles. 
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It is instructive to study the dependence of the y„A values on the minimum |At]| separation in 
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Fig. 9: (Color online) V„\ values from 0-10% central Pb-Pb collisions (points) and global fit results 
(solid lines) for 3.0 < pj < 4.0 GeV/c as a function of the minimum |At]| separation for a selection of 
Pj bins. For clarity, points are shown with statistical error bars only. For reference, a dashed line (drawn 
only in the n = 2 panel) indicates the | Atj Imin = 0.8 requirement applied throughout this analysis. 

order to observe the influence of the near-side peak. This is shown in Fig. |9] The VnA values 
rise as the pseudorapidity gap is reduced and a larger portion of the near-side peak is included 
in the correlations. At pj > 3-4 GeV/c, the peak is narrow and the curves are fairly flat at 
|At]| > 0.5. For the 3-4 GeV/c range shown in the figure, there is a discernible contribution 
from the near-side peak, but the difference does not exceed a few percent at |At] | > 0.8. 

For the first harmonic, the disagreement grows significantly as \Ar}\ is decreased, while the 
higher harmonics exhibit a much lower sensitivity to |At] |„j„,. Even if a large r] gap is applied, 
however. Fig. |9]indicates that an accurate global description of Via is still not obtained, even at 
these low to intermediate pr values. 

This behavior is representative of a general lack of consistent Vja factorization, as demonstrated 



in Fig. 10 where the global fit for n = 1 is shown for two centrality ranges. However, a rea- 



sonable fit is obtained if the pf range is restricted to smaller intervals. Fig. 1 1 shows the result 
of performing the global fit to vi (left) and V2 (right) over 0.25 < < 1.0 and 2 < pj < 4 
GeV/c separately at two different centralities. In the case of vi, a divergence occurs between 
the results obtained from the two different p^ bins, which is more prominent in mid-central 
collisions. For 0-20% V2, however, fits using the two different pf intervals lead to approxi- 
mately the same curve, supporting the observation that V2 factorizes: on average over events in 
this centrality category, a unique symmetry plane exists for the majority of all particles below 4 
GeV/c. The systematic increase of the higher pj fit compared to < 1 is likely from nonflow 
contributions on the away side, which are larger in the more peripheral centralities because of 
reduced quenching effects. Thus the observed patterns follow the expected trends with pj and 
centrality. 

The breakdown of factorization for Via does not imply that there is no real collective vi since 
the collective part may not be the dominant contribution to Via. It is therefore interesting to note 
that at low pj, the best-fit vi{GF} values become negative, as observed in hydrodynamic sim- 
ulations with fluctuating initial conditions P3ll . Although those calculations were for Au-i-Au 
collisions at 200 GeV, qualitatively similar results have been obtained at 2.76 TeV ll44]| . Esti- 
mation of the effect of momentum conservation as a correction to the coefficients prior to the 
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Fig. 10: (Color online) Global fit examples in 0-10% central (top) and 40-50% central events (bottom) 
for « = 1 . The uncertainties are represented in the same way as for Fig. [5] 
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Fig. 11: (Color online) v\{GF} (left) and V2{GF} (right) as obtained using restricted pj fit ranges 
at two different centralities. The open circles (solid squares) represent fits using only 0.25 < pj < \ 
{2 < pj < 4) GeV/c. For the more central V2 points, the two different fit ranges lead to a similar curve, 
indicating an approximate factorization. In contrast, a divergence with rising pj is observed for v\{GF}. 
For both n=\ and n = 2, the divergence is enhanced in more peripheral centrality classes. 



global fit is currently under investigation. In [|43l . such a correction amounted to a change in 
vi of about 0.01-0.02. In this analysis we have included a systematic uncertainty of 10% (p^) 
in Via to account for the bias resulting from the neglect of this correction, and the uncertainty 
is propagated to vi{GF} in the same fashion as for n > I. Further studies will be required to 
unambiguously extract the collective part of Via- 
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6 Summary 

The shape evolution of triggered pair distributions was investigated quantitatively using a dis- 
crete Fourier decomposition. In the bulk-dominated pj regime, a distinct near-side ridge and a 
doubly-peaked away-side structure are observed in the most central events, both persisting to a 
large relative pseudorapidity interval between trigger and associated particles. These features 
are represented in Fourier spectra by harmonic amplitudes, both even and odd, which are finite 
in magnitude up to approximately n = 5. These pair anisotropics are found to approximately 
factorize into single-particle harmonic coefficients for pj <A GeV/c, with the notable exception 
of Via- This factorization is consistent with expectations from collective response to anisotropic 
initial conditions, which provides a complete and self-consistent picture explaining the observed 
features without invocation of dynamical mechanisms such as Mach shock waves [|45l . 

The data also suggest that at low pt (below approximately 3 GeV/c), any contribution from the 
away-side jet is constrained to be relatively small. In contrast, for associated pj greater than 
4-6 GeV/c, the long-range correlation appears dominated by a large peak from the recoil jet. 
In this regime, when both particles are at high momenta, the anisotropy does not follow the 
;? 7 -dependent pattern followed by particle pairs at lower pj. The global fit technique provides 
a means of identifying transitions in the momentum and centrality dependence of correlations 
with respect to symmetry planes. Within the bulk-dominated region, the measurement of all 
significant harmonics provides the possibility to constrain the geometry of the fluctuating initial 
state and further understand the nuclear medium through its collective response. 
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